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We investigate the cosmological reconstruction in anisotropic universe for both homogeneous and 
inhomogeneous content of the universe. Special attention is attached to three interesting cases: 
Bianchi type-I, Bianchi type-Ill and Kantowski-Sachs models. The de Sitter, power-law and general 
exponential solutions are assumed for the scale factor in each spatial direction and the corresponding 
cosmological models are reconstructed. Moreover, for the general exponential solutions, from which 
the de Sitter and power-law solutions may be obtained, we obtain models which reproduce the early 
universe, assumed as the inflation, and the late time accelerated expanding universe. The models 
■ obtained for the late time universe are consistent with a known result in literature where a power-law 

type correction in T is added to a power-law type of f(T) for guaranteeing the avoidance of the Big 
Rip and the Big Freeze. 
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I. INTRODUCTION 
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The probably presence of an unknown form of energy in the universe, called dark energy, and confirmed by a 
£f~) • large number of observations, starting by the data of Supernovae IA in 1998 [jj, have leaded to explore the possible 
^ ' theoretical origin for this fluid. Since it is the direct responsible of the present accelerating expansion of the universe, 
a negative pressure is required which leads to a negative equation of state parameter. The most popular candidate, 
the cosmological constant, which posses a constant equation of state (EoS), p\ = — pa, can explain quite well the 
cosmological evolution. However, the open possibility that the EoS is not completely constant but evolutes dynamically 
(even crossing the phantom barrier more), and the quite large difference between the observed dark energy density 
, and the vacuum energy density predicted by quantum field theories, have leaded to explore other possibilities, as the 
existence of scalar fields, vector fields or modifications of General Relativity (GR), among others (for a review on dark 
energy candidates, see Q)- 

• • ■ In the context of modified gravities, a wide range of possibilities have been explored, being f(R) gravity probably 
the most popular one due to its simplicity since it generalizes the Hilbert-Einstein action to a more complex function 
of the Ricci scalar (for a review on f(R) gravity, see 0,l3l)- Nevertheless, other kind of theories have been suggested, 
where other curvature invariants are included as the Gauss-Bonnet gravity. In this paper, we study the so-called f(T) 
gravity, which in analogy to f(R) gravity, consists in a generalization of the action of Teleparallel gravity, a theory 
that assumes Weitzenbock connection instead of the Levi-Civita connection, which yields to a null curvature but a 
non- vanishing torsion (for a review see [5|). In this gravitational theory, the main field is represented by the so-called 
tetrads instead of the metric as in GR. This kind of theories has become very popular recently as can also explain 
the accelerated expansion of the universe with no need of dark energy, and even the inflationary epoch (see [(1- 11311) . 
Then, a wide number of aspects have been studied in the context of f(T) gravity, as its local Lorentz invariance [14j . 
static solutions , non-diagonal tetrads [l6| , or the presence of wormholes pj} , as well as other aspects 0, Qjl . Also 
a large effort has been done to study cosmological solutions for this class of theories, as well as possible cosmological 
predictions (see Refs. [20|-[23|). 

At the present work, we are interested to study some particular cosmological solutions in f(T) gravity, where the 
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appropriate action is reconstructed for each case. Specifically, the Bianchi type-I, Kantowski-Sachs (KS) and Bianchi 
type-Ill models are considered, and particularly some important solutions, such as power law and de Sitter (dS) 
expansion, or more complex ones as exponential functions for the scale factor in each direction of the space. Since 
power law and dS solutions can provide a good description for some specific phases of the universe evolution, their 
reconstruction in f(T) gravity becomes a crucial point in order to consider this class of theories as serious candidates 
for explaining the whole cosmological history. In addition, here we assume more general cosmological metrics than 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metrics, in particular anisotropic universes described by the Bianchi 
type-I, Kantowski-Sachs (KS) and Bianchi type-Ill metrics, in order to provide the most general description of the 
cosmological evolution in the context of f(T) gravity. Moreover, exponential solutions are also considered, this kind 
of expansions has become very popular recently as they may conduct the universe to a non-singular state, where some 



bounded systems may be broken. Such state suggested in Ref. [24|, and called Little Rip, has already been studied 
in /(i?) gravity (see Ref. (25|), as well as in f(T) theories (26[. Even more, the possible occurrence of a Little Rip 
has been also explored in the context of the so-called viable modified gravities (see Ref. [27]). Note that anisotropic 
cosmological metrics have been already studied in the context of GR with the presence of isotropic and anisotropic 
fluids, as well as the stability of the solutions j28l[2^|. 

Furthermore, the use of an auxiliary scalar field, in analogy to the equivalence of Brans-Dicke theories for f(R) 
gravity (see for instance Ref. [30|). is also implemented, from which may result a useful tool to reconstruct the 
appropriate action as well as for studying the properties of f(T) gravity. 

The main motivations of assuming the assumption of a model with anisotropic geometry are based on several 
physical aspects as: the famous problem of the CMB quadrupole can be solved by considering a universe with planar 
symmetry [32| where eccentricity in decoupling, generated by a uniform cosmic magnetic field whose current strength, 
B(to) ~ 10 -9 Gauss, should be close to e c iec ~ 10 -2 ; the Bianchi type models in Loop Quantum Cosmology : 33], 4 He 
abundance [34[, cosmicparallax [3^, [Hj], small anisotropic pressures (37[, cosmological solutions of the low energy 
string effective action (3a. anisotropic inflationary universe [39| and some other [Z(|. In the f(R) theory, we already 
have some good results|4ll. |42| . therefore, we propose to establish the equations here and get the first results in f(T) 
gravity, for the Bianchi type-I, type-Ill and KS models. 

Then, the paper is organized as follows: in sectionflll the basic concepts of /(T) gravity are introduced. In section 
IIIII , the equations for general Bianchi type-I, type-Ill and Kantowski-Sachs (KS) models are deduced in a particular 
coordinate system and diagonal tetrads. Section HVl deals with the reconstruction of the f(T) action for some relevant 
solutions, and where several techniques are considered, including a kind of scalar-tensor theory for torsion gravity. 
Finally, section [V] is devoted to the conclusions and discussions on the results found in the paper. 

II. PRELIMINARY DEFINITIONS AND EQUATIONS OF MOTION 

As previously mentioned, the f{T) theory of gravity is defined in the Weitzenbock's space time in which the line 
element is described by 

dS 2 = g^dx^dx" , (I) 

where g^ are the components of the metric which is symmetric and possesses 10 degrees of freedom. One can describe 
the theory in the spacetime or in the tangent space, which allows us to rewrite the line element (fT|) as follows 

dS 2 = g^dx^dx" = r] ij 6 i e j , (2) 
dx" = e i ^6 i , 6 i = e\dx» , (3) 

where rjij = diag[\, —1, —1, —1] and e^e 1 v — 8^ or e i „ — 5\. The square root of the metric determinant is given 
by ^/~— g = dot [e l ] = e and the matrix e a M are called tetrads and represent the dynamic fields of the theory. 
By using theses fields, one can define the Weitzenbock's connection as 

T^ u = e i a a v e i ^-e\d„e t a . (4) 

The main geometrical objects of the spacetime are constructed from this connection. The components of the tensor 
torsion are defined by the antisymmetric part of this connection 

T a yLv — — r"„ = ej" (p^e 1 v — d v e l ^ . (5) 

The components of the contorsion are defined as 

K^a = -\{T^ a -T^ a -T a n . (6) 
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In order to make more clear the definition of the scalar equivalent to the curvature scalar of RG, we first define a new 
tensor , constructed from the components of the tensors torsion and contorsion as 

ScT = \ (k^* + ^% - ^ T %) ■ (7) 

We can now define the torsion scalar by the following contraction 

T = T^ v Sr . (8) 
The action of the theory is defined by generalizing the Teleparallel theory, as 

S = [ e [/(T) + Cuatter] dtx , (9) 



where f(T) is an algebraic function of the torsion scalar T. Making the functional variation of the action ([9]) with 
respect to the tetrads, we get the following equations of motion [lil TlH [2lj 

S^ p d P Tf TT + [e-%d p (ee, a SJ") + T* X ^SJ X ] f T + U^f = 4tt7^ , (10) 

where 7^ is the energy momentum tensor, f T = df(T)/dT and f TT = d 2 f(T)/dT 2 . By setting f(T) = a x T + a , 
the equations of motion (fTU| are the same as that of the Teleparallel theory with a cosmological constant, and this is 
dynamically equivalent to the GR. These equations clearly depend on the choice made for the set of tetrads . 

The contribution of the interaction with the matter fields is given by the energy momentum tensor which, is this 
case, is defined as 

7/ = diag(l,-oj x ,-uj y ,-uj z ) p , (11) 
where the u>i (i — x, y, z) are the parameters of equations of state related to the pressures p x , p y and p z . 

III. FIELD EQUATIONS FOR BIANCHI TYPE-I, TYPE-III AND KANTOWSKI-SACHS MODELS 

Let us first establish the equations of motion of a set of diagonal tetrads using the Cartesian coordinate metric, 
for describing models of Bianchi type-I, type-Ill and Kantowski-Sachs (KS). We propose to start with the Bianchi 
type-Ill case, from which Bianchi type-I and KS can be recovered. For the Bianchi type-Ill case, the metric reads 

dS 2 = dt 2 - A 2 (t)dx 2 - e- 2ax B 2 (t)dy 2 - C 2 {t)dz 2 , (12) 

where a is a constant parameter. Note that the Bianchi type-I is recovered by setting a — from the Bianchi type-Ill, 
while KS is recovered when one takes a = and B(t) = C{t). Let us choose the following set of diagonal tetrads 
related to the metric (|T2"|) 

[e a J = diag [l, A, e~ ax B, C] . (13) 

The determinant of the matrix flT3l is e = e~ ax ABC . The components of the tensor torsion for the tetrads (TlB"]) 
are given by 

T 1 — — T 2 - — T 2 -a T 3 -— (U) 
01 ~ A 02 ~ B 21 ~ 03 _ C ' 

and the components of the corresponding tensor contorsion are 

K Q \ = \ , K°\ = | , K^ 2 2 = £ , ^°3 3 = g . (15) 
The components of the tensor S^, in 0, are given by 
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By using the components (IT4l) and (IT6|) . the torsion scalar ([8]) is given by 

T = -2 

The equations of motion are given by 
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In the particular case where f(T) = T — 2A, the equations (fT8")) - (l2U)) are identical to that of the GR 43]. The equation 
of constraint (|23|) appears in both the GR as in f(R) gravity [42j. But here we have a second equation of constraint 
(|22[) . which appears as a generalization of the previous one, because here we have a contribution of a term of second 
derivative of the function f(T) with respect to T. 

By setting a — 0, the Bianchi type-I case is recovered and the equations of motions read 
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The equations of motion corresponding to KS model are obtained by setting a = and B = C ', yielding 
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ltorpy = f + 2f T 
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In the next section we will perform the reconstruction scheme of the action of the system for some particular cases. 



IV. RECONSTRUCTING f(T) GRAVITY IN INHOMOGENEOUS UNIVERSES 

Let us now consider the reconstruction of the f(T) action for some particular solutions of the class of metrics 
explored in the previous section. Specifically, we consider solutions of the type of de Sitter, power law evolutions and 
exponential solutions. Note that de Sitter and power law solutions have been widely explored in other contexts of 
modified gravity, as f(R) and Gauss-Bonnet gravities (see Ref. [44|), since they can provide a well description of the 
cosmological evolution along its particular phases. 

Let's start by considering for simplicity Bianchi type-I and Kantowski-Sachs (a — 0) metrics. Then, the conservation 
equation for the energy momentum tensor (111)) can be easily obtained, 

p + (H x + H y +H z )p + H xPx + H yPy + H z p , = , (31) 

where we have defined H x = ^ H y = -g H z — ^. We can now analyse de Sitter, power law solutions and exponential 
expansion in Bianchi type-I metric by one side, and Kantowski-Sachs metric by the other, where B = C that implies 

Py = Pz- 



A. De Sitter solutions 



De Sitter solutions are well known in the context of cosmology since the current epoch, where the universe expansion 
is being accelerated, can be described approximately with a de Sitter solution. This kind of solutions consists on an 
exponential expansion of the scale factor, which yields a constant Hubble parameter. In the case of Bianchi type-I 
and Kantowski-Sachs metrics (a = 0) in (|12|). we may assume an exponential expansion for each spatial direction, 

A = A a c at B = B a c bt C = C e ct , (32) 

and the rates of the expansion for each direction can be defined as, 

ABC 
Hx = = H x0 H y = — = Hy H c = — = H z0 , (33) 

where {H x q — a, H y o — 6, H z q = c} are constants. The torsion scalar defined in (1171) is given by, 

To = —2 (H x oH y o + H x0 H z0 + H y0 H z0 ) . (34) 

Then, by assuming p x = p„ = p z = p and an equation of state p — wp, the conservation equation (1311) can be easily 
solved for the ansatz (1321). 



p = p e-^o+H y0 +H M )(l+ W )t (35) 

Hence, the field equations ([M|) -([2T | become, 

16 7 r (0o e- ( ^ 0+ ^ 0+ ^ o)(1+lu)t = /(To) + 4 [H x0 H y0 + H z0 (H x0 + H y0 )} f T (T ) , (36) 

-16irwpoe-( H * 0+H y° +H ^ 1+w '> t = f(T ) + 2(H y0 + H z0 )(H x0 + H y0 + H zQ )f T (T ) , (37) 

~l&i:wp (s e^ H ^ +H y a+H ^ 1+ ^ t = /(T ) + 2(H x0 + H z0 )(H x0 + H y0 + H z0 )f T (T Q ) , (38) 

-167Twp Q e-^ 0+H y" +H ^ ( - 1+w '> t = /(T ) + 2(H x0 + H y0 )(H x0 + H y0 + H z0 )f T (T ) . (39) 

Note that the only possible solution in the presence of a perfect fluid is one with w = — 1 as the r.h.s. of equations 
(|3"6")) - (|39p is independent of time, according to the expression of the scalar torsion for a pure de Sitter solution (|34l) . 
unless H x q + H y o + H z q = 0, which would imply a decelerating expansion in a particular direction, being < 0. 
Moreover, for a particular f(T) action, the system of equations (|3~4")) - (j3l?]) reduces to an algebraic system of equations 
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for the variables {H x o, H v q, H z q}. Since the system of equations (|34]) - (j39|) are composed by four equations, while there 
are only three variables, the above 4-equations system has to be reduced. However, even in the case of Kantowski- 
Sachs metric, where B(t) = C(t) —> H y0 = H z0 , the system (|34"1) - (|39|) still posses three independent equations with 
two variables. Hence, the only possible solution imposes, 

A(t) = B(t) = C(t) ->■ H x0 = H y0 = H z0 = H Q , (40) 

And the metric ()12j) reduces to the well known Friedmann-Lemaitre-Robertson- Walker metric with an exponential 
expansion, A(t) = A$e Ho *. Hence, the only solution for a pure de Sitter expansion in Bianchi type-I and Kantowski- 
Sachs metrics gives a FLRW universe |4&J, and the system of equations (j3"4")) - (l3l?l) reduces now to a unique independent 
equation, 

IGttpo = /(T ) + \2Hlf T {T ) . (41) 

Then, the roots of the algebraic equation (j4Tj) give the de Sitter points of a particular f(T) action. In order to 
illustrate such possibility, let us consider the action, 

f(T) = (-T)" , (42) 

where n is a real constant. Then, the equation (|41j) is rewritten as, 

16^ = (1 - 2n)(6ff 2 )" , (43) 

whose solution is given by, 

r2 1/167TAA 1 /" 
6 \ l-2n) 



H^ = -A ■ (44) 



Hence, the only physical solution (p ,Hq > 0) imposes n < 1/2. Then, the de Sitter solution is a direct consequence 
of the energy density po, which can be interpreted as a cosmological constant according to the condition imposed 
above for its equation of state, w = — 1. Nevertheless, in vacuum the equation (l43t reduces to = (1 — 2n)(QH^ ) n , 
whose only solution is given by n = 1/2, rising to f(T) = y—T that posses an infinite number of de Sitter points. 
Moreover, we may consider in vacuum the action, 

/(T) = CxT + C 2 (-T) n , (45) 

where {C\, C2} are the coupling constants and n is a real constant. The field equation (1411) in vacuum yields, 

= d&Hl +C 2 (1 - 2n)(6# 2 ) n . (46) 



So the roots of this equation give the dS points allowed by the class of theories expressed in (|45|) . Note that now, 
the exponential expansion is a direct consequence of the action instead of the contribution of a kind of cosmological 
constant as in the case shown above. For instance, n = 2, it yields the solution, 



While for higher powers of n, more de Sitter points can be obtained for the action (|4"5"j) . Note that in f(R) theories, 
dS points constitutes the critical points of the dynamical system, which may be (un)stable, and could explain both 
the inflationary and dark energy epochs (see [4g]), which may be the case also in f(T) gravity. 



B. Power law solutions 



Let us now explore a cosmological evolution described by a power law in each direction of the space expansion. In 
such case, the scale parameters for the Bianchi type-I and Kantowski-Sachs metric (p~2|) . where we set (a — 0), can 
be expressed as, 

A(t) = A Q t a , B(t) = B t b , C(t) = C t c , (48) 
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where {a,b,c} and {Aq, Bq,Co} are constants to be determined by the field equations, and initial conditions respec- 
tively. The expansion rates are given by, 

a b c , . 

H x = - , H y = - , H z = - . (49) 

While the expression for the torsion scalar (fTT)) yields, 

_ „ / ab ac bc\ ,„„ s 

T = 2 {¥ + ¥ + *)- (50) 

Then, introducing the above quantities in the field equations (|M)) - (I?7|) , we get the following system of differential 
equations in f(T), 

16ttp(T) = f(T) - 2T/ T (T) , (51) 

-W«MT) - f(T) + {b+ fl~ a ~ b ~ C) TMT) + 2 . j 6 + C l , T 2 fxT(T) , (52) 

oc + a(o + c) oc + a(o + c) 

-lft™,p(T) = /(T) + {a+ ? { l-°- b ~ C) TMT) + 2 j a + C > T*f TT (T) , (53) 

oc + a(o + c) oc + a(o + c) 

-16™ P (T) = /(T) + (0+ . h) ^" fl " h " C) r/ r (r) + 2 ^ | a + 5 | - TVrr(T) , (54) 

oc + a(o + c) be + a(b + c) 

where we have assumed for simplicity that p x = p y = p z = p and an EoS p — wp. Hence, the system (|5 1 [) - (|54[) is a 
set of differential equations in f(T) with the torsion scalar T as the independent variable. 

Firstly let us consider the vacuum case, or in other words, the homogeneous part of the first equation (|51[) . which 
becomes f(T) — 2T/y(T) = and whose solution is given by, 

f(T) - dV=T , (55) 

where C% is an integration constant. In order to satisfy the rest of the equations ([52^ -([54 |) . the condition a = b = c must 
be imposed, so that the Hubble parameters (|49|) reduce to the usual FLRW cosmology reproducing power law solution. 



In the presence of a kind of isotropic perfect fluid p = wp, we can first solve the continuity equation pip in order 
to obtain p = p(T), which yields, 

(q + 6+c)(l+ni) 

P = P ot- {a+b+c ^ = Po [- -, T ) ' . (56) 

\ 2{ab + ac + be) ) 

Hence, the solution for the set of equations ([ST|) - ([M)) is given by f(T) = fh(T) + f p (T), where fh(T) is the solution 
of the homogeneous equation, which coincides with the vacuum solution (I55p . and f p (T) is the particular solution. 
Then, by using (|56|) in the equation (|5Tj) . the particular solution can be easily found, 

f p (T)= X T ^ l , (57) 

where x is a constant given by 

2 i-(l+w)(a+b+c)/2 7TpQ 



[-1 + (1 + to) (a + b + c)] [-be - a(b + c)] <1+ ™ ) 



(58) 



Note that the condition (1 + to) (a + b + c) = 2 n with n being any natural number, has to be imposed in order to 
avoid a complex gravitational action that would lacks any physical sense, recall that T < for an expanding universe 
according to (|50j) . In order to satisfy the complete set of equations (|51"1) - (|5~4"|) . we introduce the solution ([57]) into the 
field equations (f52"]) - (p>4")) , and the following solutions for the parameters {a, b, c} are found, 

i. c = 1 ~ wt ^ +b ^ > , where to 7^ 0. This provides an anisotropic solution in f(T) gravity with A(t), B(t) and C(t) 
being different functions in (|4"5|) . and recalling that the perfect fluid assumed is an isotropic fluid. This does not 
hold in GR or Teleparallel Theory (TT) but is possible here due to the presence of second derivatives of the 
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function f(T) with respect to the torsion scalar T in ([5Tj) - ([54f . Note that field equations may be rewritten as 
the usual equations in TT, 

H x H y + H X H Z + H y H z = 16tt(p + p f{T) ) , -H y - H 2 — H z — H 2 Z -H y H z = 8tt(w P + p x f(T) ) , 



-H x - H 2 X -H z - H 2 Z -H X H Z = 8tt(wp + p y f{T) ) , -H x - H 2 X - H y - H 2 - H x H y = 8tt(w P + p z f(T) ) . (59) 

Here, the extra terms coming from /(T) are defined as an energy density P/(t) and pressures {P/(t), -P/(t)' P/(t)}' 
which are the origin of the anisotropic evolution. In this case, we have to fix C\ — in (|55[) in order to satisfy 
the whole system. 

ii. a = b = c. The cosmological evolution expressed by [48] reduces to a FLRW metric as in the homogeneous part 
of the equations, so that C\ ^ 0. 

iii. c = —ab/(a + b). In spite of this satisfies equations (|51"1) - (|M)) once (|57|) is substituted in the equations, it gives 
T = 0, and the r.h.s of equations (|51"1) - (|M)) become null while the l.h.s. are not, since p = p(t) as given in (|56l) . 
so this is not a real solution. 



Therefore, we have obtained a complete set of power law solutions for Bianchi type-I universe and Kantowski-Sachs 
metrics in the context of f(T) gravity. Nevertheless, the action is clearly dependent on the EoS parameter w. Note 
also that in vacuum, the only possible solution reduces to a FLRW metric. 



C. General exponential solutions 



In this subsection we consider a more general exponential expansion for each spatial direction by 

A = A eM t] , B = B e 9 * (t) , C = C e s * (t) , 
where the function gi(t) is assumed as 

9i(t) = hi(t)\n(t) , i = x,y,z , 



(60) 



(61) 



and Aqi Bq and Co are positive constants. Note that the previous cases, the de Sitter solutions and power law solutions 
can be recovered from this one by setting hi(t) = ait/(hx(t)) and hi{t) = aj, respectively, where {a^} = {a, b, c}. In 
what follows, we will use an adiabatic approximation for the expansion in each spatial direction and neglect the 
derivatives of hi(t), i.e., setting (hi ~ hi ~ 0). The expansion rates in this case are given by 



H r 



Thus, the torsion scalar (fTTf becomes 



h x {t) 



H y = 



hy(t) 



H z = 



h z (t) 



h X {t)hy{t) h X {t)h Z {t) hy{t)h Z {t) 



t 2 



f 2 



The acceleration in each direction is given by 



A= K{K ~ l) A, B = 



_ K{hy - 1 
t 2 " ' " t 2 " ' " t 2 



B, C= h '« l '- 1) C 



(62) 



(63) 



(64) 



Since A, B and C are positives, the acceleration is guaranteed in each direction when hi > 1, while for < hi < 1, 
the universe is in deceleration. 
The simplest example of hi (t) is 



hi(t) = 



Hin ~r ' H out 
l + qt 2 



qt 2 



(65) 
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where ha n , hi tOUt and q are positive constants, and the q is assumed to be enough small in order to make hi{t) varying 
slowly. Thus, the torsion scalar is always negative. From (|65[) . we see that at early time t = 0, hi — > h iin and for late 
universe, hi —> h iout . By using (|65|) and ff63]) . one gets the following equation 



g 2 Tt 6 + (2qT + 2q 2 X) t + (T + 2qY) t 2 + 2Z = , 

h x ou i hy out ~\~ h x outhz out ou f h z ou t , Y — h x i n hy ou t -\~ h x out^y in ^x in ^z out H~ 

h x out^z in H~ hy in^z out H~ ^y out^z in j ^ hx in^y in ^xin^zin H~ ^yin^zin 

whose solutions read 

i 2 = {*„(T), *±(T)} , 
* (T) = a 1/3 +^ a- 1/3 + /3i *±(T) = e T2 " /3 a 1/3 + e ±2i7r/3 a~ 1/3 + ft , 

■ /HT 9 T(8X - 6Y) + 4g 2 X 2 + T 2 „ 9nX + 9T 

«2 , PO = » 9rTl9 , Pi = - 



(66) 



(67) 



g 2 T 2 
1 



a 2 



9g 2 T 2 

27qT 2 Z 2 + [{-36q 2 TX - 36qT 2 )Y + 16g 3 X 3 

4g 3 X 2 - 8g 2 TX + 8qT 2 )Y 2 + (-Aq 2 TX 2 - 8gT 2 X + 2T 3 )F - qT 2 X 2 - 2T 3 X 



21qT 
8q 2 TY 3 

qT 2 (27 Z - 18Y + 15X) + q 2 T(24X 2 - 18XY) + 8q 3 X 3 - T 3 



3qT 

2 TX 2 + 30qT 2 X ~ 2T 3 ]Z 



27q 3 T 3 

We see from (|67|1 that there are one real positive solution, ^o(7 1 ), and two complex solutions ^>±(T). By using 
([62"T) - (|64[) the system of equations of motion (|24 |) -([27 |) becomes 



16* p = / - 2Tf T , 
-16np x = f + 2f T 

-IQirpy = f + 2f T 



-16tt Pz = f + 2f T 



{hy ~ 


Yh z f 


- (hy + h Z ) 


- h y h z 






*o(r) 




(h x - 


Vh z ) 2 


- (h x + kg) 


- h x h z 






*o(T) 




(h x - 




- (h x + h y ) 


h X hy 



T 
2 
T 
2 



*o(T) 



h v + hz 

*o(T) 
h x -f- /iy 



TJtt 
TJtt 



hi = 



h iout q^o(T) 



l + q* (T) 



(68) 
(69) 

(70) 

(71) 

(72) 



i. A special case 



Here, we assume that the expansion rates are equal in the three spatial directions (h x — h y = h z ), and the system 
([7T|) reduces to 



16* p = / - 2Tf T , 
-I6np x = f~2f T (t 

Px = Py = Pz , 



' -2T 
3*o(T) 



AT* 



I -2T 
3*o(T) 



TT 



(73) 
(74) 



which means that the assumption of having the same rate in the three spatial direction leads to an isotropic matter 
content. By combining (|73]l and ([74]). one gets 



4T 4 



/ -2T 
3*o (T) 



/tt — 2/t 



(1+cj)T + 



' -2T 

3* (r) 



+ (l + w)/ = 



(75) 



where we used the barotropic equation p x = up. Let us consider an asymptotic analysis, looking for the early universe 
(small time) and the late time one (large time), for which the function hi(t) yields ha n and hi OUtl respectively. Thus, 
for the early universe Eq. (|75p reduces to 



4T 2 / TT + 6T/ T [3h im (tu + 1) - 1] - 3h iin (oj + 1)/ = , 



(76) 



h% in — h x i n — hy — h z i n , 
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whose general solution reads 



f(T) = C a T+ + C^- , 7± - 5 - gWl + ") ± V 25 - TgWl + gj + 81/& n (l + ^) 2 ^ (??) 

where C3 and C4 are integration constants. From (|77l) . by writing the radicand as [5 — 9h X i n (l + uj)] 2 + 12h X i n (l+u), 
one observes that, for any uj > — 1, 7+ > and 7_ < 0. Moreover, in this context of asymptotic analysis, we observe 
from (|63p that for small t, the torsion scalar T is large, while for large t, the torsion is small. Thus, for small t with 
hxin > 1) corresponding to the inflation, the algebraic expression of f(T) is given by 

/(T)=C 3 TT+ . (78) 

Since h x (t) reduces to h xou t in the late universe, the model corresponding to the late accelerated universe can be 
obtained by replacing h x i n by h xout . Precisely, for large t, the torsion scalar is small, and for h xou t > 1, the dominate 
term in (|77l) . corresponding to the model of late time universe, is 



f(T , rT1 > , 5 - 9/w(l +u)± y/25 - 7&h xout (l + + %lhl out {\ + ujY 

Ji 1 ) = , 7± = ^ ■ (< 9 ) 

This model is equivalent to the teleparallel gravity for C4 = 1 and h xout = 2/(5 + 5w). It is easy to see from this that, 
for any ordinary matter, i.e., uj > 0, one gets h xout < 1, meaning that the teleparallel gravity without cosmological 
constant cannot provide the late acceleration of the universe (remembering that the acceleration is guaranteed for 
h xou t > 1, and < h xout < 1 characterising a decelerated expanding universe). Thus, the contribution of the f(T) 
terms plays the role of the dark energy. 

Looking for the expression of f(T) for large cosmic time t, i.e., the expression (|77|) (replacing h xin by h xout ), a 
similarity can be observed with a result of Bamba et al in |26j. In this work, they undertook f(T) theory in the 
FLRW metric, first assuming a power-law expression for f(T) in the Eq. (4.6) for investigating what type of finite 
future time singularities may appear. Later, they introduced a correction term, still in the form of power-law, their 
Eq. (4.22), in order to analyse the possible avoidance of the singularities. Then, they obtained the global expression 
(4.23) of [26j]. Note that this expression is equivalent to our Eq. (|77l) . Moreover, they shown in their "TABLE II" 
that the Big Rip and the Big Freeze can be removed if the product of the exponents is negative. This is exactly our 
result, since 7^ and are always positive and negative, respectively ( 7 ^_ are obtained from j± by replacing h x in by 
h xo ut)- This means that for both Bianchi type-I and KS, if the expansion of the universe occurs with the same rate 
in all direction, models that can realize the late time accelerated expansion of the universe, and able to prevent the 
Big Rip and the Big Freeze can be reconstructed. 



2. Using an auxiliary scalar field 

In this subsection, we would like to use the reconstruction scheme, for which an auxiliary scalar field is introduced 
in the action of the theory. By this way, the functional form of f(T) can be found through two other scalar functions 



P and Q (for more clarification on the method, see the case of f(R) gravity [46j and other of f(T) gravity [16|). Here, 
we return to the general system (p5|) - ([72|) . where we introduce the auxiliary scalar field 4>, by defining the algebraic 
function f(T) by 

f{T) — P{4>)T + Q{4>) . (80) 

By using (15(H) and varying the action © with respect to the scalar field <fi, one gets 

P'{cj>)T + Q'{cj>)=Q , (81) 

which may be solved with respect to (f>, yielding cf> — <f>(T). Here, the prime (') denotes the derivative with respect to 
4>. By using (|50"|) . one obtains 

f T (T) = P(0(T)) , Jtt(T) = PrWH) ■ (82) 
Making use of ((821) . one can rewrite the system ([68|) -([7T |) as 

16^ = P{0)T + Q(<f>)-2TP{0) , (83) 

-16^ = PW + QM+2PW [ (^±M!^+M^Mi _ Z\ 

k ^f)TP T (*), (84) 
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Wwpy = P(4>)T + Q(^)+2P(<p) 
hx + h. 



(h x 



(K 



h x h z 



P(4>)T + Q(4>) + 2P(4>) 



(h x + h y ) 2 - (h x + h y ) - h x h y 



*o(T) 



*o(T) 

Let us consider that p y = p z . Then, by equating 

TP T ((f>) = - 
Using (J87J), the system (|83 ]) -([86 ]) reduces to 

167T/5 

-Idirpx 

Px 



r tPtW 



with 



one gets 



1) 



Q(4>) - TP(4>) , 
Q(4>)+TP(4>) , 

Py = P Z , U X = LOy 



By using (|88|) . one can determine Q((f>) as 



which, substituted in ([M]) and yields 



Q(cj)) = 16tt p+TP((j>) 



1Q-k(w x + l)p = 2TP(4>) 



(85) 



(86) 



(87) 



(88) 
(89) 
(90) 



(91) 



(92) 



As one can redefine the scalar field properly, we may choose 4> = t, then P(4>) = P(t) = P(T). Note that in the case 
where lu x — — 1, one gets P(T) = 0, and the algebraic function f(T) = Q(i) = 16irp = Q(T). But this requires to 
have the complete expression of p depending on T, i.e., solving the equation of continuity. Let us consider in general 
the case where uj x ^ —1 and try to solve the equation of continuity, which will help us to determine P(T) and Q(T), 
leading to the reconstruction of the algebraic function f(T). We choose the case of KS geometry (h y = h z ), where 
the equation of continuity is written as 



p + p(l+cu x )(H x + 2H y ) =0 



which can be solved, giving 

Pit) 
Gi(t) 



C 5 exp[G 1 (t) + G 2 (i)] , 

(1 + bJ x )(h x i n + 2hyi n ) 



G 2 {t) = s/q[(l + u x )(h x 



2h 



y xn 



2h vout )] arctan( v / gt) 



(93) 

(94) 

(95) 
(96) 



where C5 is an integration constant. Making use of (|67]). one can cast p(t) in terms of T. Also, from Eq. (|92l) . one 
obtains P(t) in functions of T, and from (|9Tj). one gets Q(t) in terms of T. Therefore, we can reconstruct the algebraic 
function f(T), given in (j80| . as 



f(T) 



-16irC3Ui x exp 



,9i 



gi arctan 



9*ocn) 



(97) 



2h„ 



2h 



y out 



where gi = (1 +uj x ){h xm + 2h yin ) and g 2 = ^/q(l +u x )(h xin - h xo 

In principle, with (|97|) . some cosmological models can now be reproduced. Let us focus our attention to the early 
universe (which may be characterized by the inflation), and the late time universe (characterized by its accelerated 
expansion) in the three spatial direction. 

At early time, i.e., for small t (or large T), ^o(T) is very small, and the corresponding algebraic function is 



f(T) = -16itC 3 uj x exp 



.91 



V*o(T) 



(98) 
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At late time, the time t is large (corresponding to small torsion scalar and large ^o(T)), and the algebraic function 
reads 



f(T) = -16ttC 3 w x exp 



g 2 arctan (y^qW^T) 



(99) 



Thus, we see that models corresponding to the inflation and the late time accelerated universe can be reconstructed 
within KS metrics where the matter content is partially isotropic (p y = p z ). Since in this work we fall in the situation 
where oj x = u) y = us y one could use the WMAP data and try to check if they fit with this anisotropic model. Because 
of having ultimately p x = p y = p z and h y = h z , one can just use the first two equations of motion of KS case, i.e., 



In order to put out the contribution of the modified part of the algebraic function f(T), we cast it into 
the form f(T) = T + j(T), as the teleraparallel term plus the modified additive algebraic function j(T). Thus, the 
equations (|2"51) and become 



(100) 
(101) 



8-rrpeff = 

SlTPxeff = 

where p e ff and Pxeff are the effective energy density and effective pressure, respectively, and defined by 

2 




Peff = P 



j + 4ja 



1 

16tt 

Pxeff = Px + { j + 4/7 




4—fjTT 



By dividing f|101[> by (|100|) . using w eff =p xef f/p eff) one gets 



2 - ih y 
h„ + 2h q 



(102) 
(103) 

(104) 



Since the observational data from the WMAP project do not depend on the spatial direction, i.e., are not based on an 
anisotropic geometry, we have to impose h x = h y in (I104[) . for getting a suitable cosmological feature. Thus, assuming 
that the accelerated expansion of the universe is realized with the same rate in all directions, h x = h y = h z > 1, and 
that the universe is essentially filled by the dark energy, where we can neglect the usual matter content such that 
Lo e ff ~ i^DE- Therefore, using (|104p . we have the possibility to regain the well known range of values allowed by the 
current 7- year WMAP data for the parameter of the equation of state of the dark energy, u>de = — 1.1 ± 0.14 WMAP 

B 



D. On Bianchi type-Ill solutions 

In this section, we propose to present some comments on Bianchi type-Ill solutions. This case is quite exceptional 
since we do not have the freedom of making cosmological reconstruction as in the case of Bianchi type-I and KS, due 
to the constraints equations (|22[) and ([23| . 

From (|23|) . since the parameter is different from zero and the algebraic function cannot be a constant, one gets 

i - ! • < 105 » 

which, injected in (|22|) leads to 

Tf TT = , (106) 



meaning that one has T = or Jtt — 0. The first case, T = implies that one has a constant torsion scalar, i.e., 
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where K is positive constant. Let us consider A — C n , with n bigger than zero or less than —2. Thus, Eq. (|107|) can 
be solved yielding 



C(t) = Cpexp ± 1 



K 



n(n + 2) 



leading to 



A(t) = B(t) = (C£)"exp ±n x 



K 



n(n + 2) 



(108) 



(109) 



where C' is a positive constant. It is important to note that in order to guarantee the expansion of the universe, one 
need to have 



A(t) = B(t) = 



(Co)"exp(n v /^ 



for n < —2 
for n > . 



(110) 



In this case, we see that the rate of expansion is constant for the three spatial direction: this is the de Sitter 
universe. 

Now we can perform the reconstruction of the algebraic function f(T). One can cast Eqs. <|18j > -(|2i p in the following 
system 



Wnp = f + 4f T [K 



-16ttuj xP = f + 2Kf T 



-16iroj z p = f + AKf T 



—) 

2 A 2 J ' 

2n 2 + 3n + 1 
n(n + 2) 
(2n+l\ n a 
\n + 2 



- 2 — f T 
A 2 



Px = Py ■ 

By combining (|111[) with (|113[) . one can eliminate the term containing a, obtaining 

-16^ + 1)^=^11/,. 

n + 2 



(111) 

(112) 

(113) 
(114) 

(115) 



The energy density p can be eliminated by combining (I112p with (|115p yielding the following differential equation 

2K [2n(n- 1)oj x - (u z + l)(2n 2 + 3n+ 1)] f T -n(n + 2){oj z + 1)/ = , (116) 
whose general solution is 

f(T) = C 6 exp[R(n)T] , (117) 

R(n) = n(n + 2) (i^ + l) 

1 ' 2K [2n(n - 1)oj x - (uj z + l){2n 2 + 3n + 1)] ' 

where Cg is an integration constant. Note that for n = 1 and oj x = uj z1 Eq. (|40[) is recovered. 

The second case from (|106[) . /tt = 0, implies that fx is constant, that we choose to be minus two times the 
cosmological constant A, then, f(T) is written as 



f(T) = T — 2A 
which is the teleparallel gravity with cosmological constant. 



(118) 



V. CONCLUSION 



Along the paper, the Bianchi type-I, Kantowski-Sachs, and Bianchi type-Ill metrics have been studied in the context 
of f(T) gravities. Particularly, we have shown the reconstruction of some important cosmological solutions, obtaining 
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the corresponding f(T) action. We have assumed initially a particular choice of coordinates and tetrads, specifically 
cartesian coordinates and a diagonal set of tetrads have been imposed in order to avoid the well known constraint 
f TT = 0, which reduces trivially to the action of teleparallel gravity (see Ref. p7j). 

Then, several important cosmological solutions have been considered. In particular, dS solutions, where the scale 
factor is an exponential function of the cosmic time, has been considered for Bianchi type-I and Kantowski-Sachs 
metrics by imposing a particular exponential expansion in each direction of the space. We have shown that the only 
possible solution turns out to the FLRW metric, such that no possible dS anisotropic evolution can be found in f(T), 
unless one considers an anisotropic fluid. Nevertheless, in the case of power law solutions, we have found that in the 
presence of a perfect isotropic fluid, an anisotropic cosmological evolution can be found for a particular choice of the 
action f(T), while in vacuum the action reduces to FLRW metric. 

Moreover, we have extended the cosmological reconstruction scheme to a general exponential solutions, from which 
the above de Sitter law and power law solutions are particular cases. We have assumed an adiabatic approximation for 
the expansion in each spatial direction. We undertook two cases, a special case and a second where an auxiliary field 
is used. In the both cases, we shown that the models can realize the early accelerated universe, characterized by the 
inflation, and the late time acceleration of our current universe. In the special case, the model presents an interesting 
aspect because it ensures the avoidance of the Big Rip and the Big Freeze. In the case where the auxiliary field is 
used, the model corresponding to the late time accelerated universe fits with the 7-year WMAP data, confirming the 
consistency of the result . 

The Bianchi type-Ill case presents some constraints from which only two forms of the algebraic function f(T) can 
be obtained. The first is the well known teleparallel gravity with cosmological constant, and the second is a de Sitter 
type solution. 
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